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METHOD AND APPARATUS FOR AUTOMATIC 
FREQUENCY CONTROL IN A CDMA RECEIVER 

BACKGROUND 

5 This invention relates generally to a method and system for controlling the reference 

frequency in a radio receiver. More particularly, this invention relates to a method and system 
for estimating a frequency offset between a carrier frequency of a transmitter and a local 
reference frequency of a receiver in a communication system. 

FIG. 1 A is a block diagram of an exemplary cellular radiotelephone system, including an 

10 exemplary base station (BS) 110 and a mobile station (MS) 120. Although denoted a "mobile 
station", the station 120 may also be another type of remote station, e.g., a fixed cellular station. 
The BS includes a control and processing unit 130 which is connected to a mobile switching 
center (MSG) 140 which in tum is connected to a public switched telephone network (PSTN) 
(not shown). General aspects of such cellular radiotelephone systems are known in the art. The 

15 BS 1 10 handles a plurality of voice channels through a voice channel transceiver 150, which is 
controlled by the control and processing unit 130. Also, each BS includes a control channel 
transceiver 160, which may be capable of handling more than one control channel. The control 
channel transceiver 160 is controlled by the control and processing unit 130. The control 
channel transceiver 160 broadcasts control information over the control channel of the BS or cell 

20 to mobiles locked to that control channel. It will be understood that the transceivers 150 and 
160 can be implemented as a single device, like the voice and control transceiver 170, for use 
with control and traffic channels that share the same radio carrier. 

The MS 120 receives the information broadcast on a control channel at its voice and 
control channel transceiver 170. Then, the processing unit 180 evaluates the received control 

25 channel information, which includes the characteristics of cells that are candidates for the MS to 
lock on to, and determines on which cell the MS should lock. Advantageously, the received 
control channel information not only includes absolute information conceming the cell with 
which it is associated, but also contains relative information conceming other cells proximate to 
the cell with which the control channel is associated, as described for example in U.S. Patent 

30 No. 5,353,332 to Raith et al, entitled "Method and Apparatus for Communication Control in a 
Radiotelephone System". 

Modem communication systems, such as a cellular radiotelephone system of the type 
described above and satelhte radio systems, employ various modes of operation (analog, digital, 
dual mode, etc.) and access techniques such as frequency division multiple access (FDMA), 
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time division multiple access (TDMA), code division multiple access (CDMA), and hybrids of 
these techniques. 

In North America, a digital cellular radiotelephone system using TDMA is called the 
Digital Advanced Mobile Phone System (D-AMPS), some of the characteristics of which are 
5 specified in the TIA/EIA/IS-136 standard published by the Telecommunications Industry 

Association and Electronic Industries Association (TIA/EIA). Another digital communication 
system using direct sequence CDMA is specified by the TIA/EIA/IS-95 standard. There are 
also frequency hopping TDMA and CDMA communication systems, one of w^hich is specified 
by the EIA SP 3389 standard (PCS 1900). The PCS 1900 standard is an implementation of the 
10 GSM system, which is common outside North America, that has been introduced for personal 
communication services (PCS) systems. 

Several proposals for the next generation of digital cellular communication systems are 
currently under discussion in various standards setting organizations, which include the 
International Telecommunications Union (ITU), the European Telecommunications Standards 
15 Institute (ETSI), and Japan^s Association of Radio Industries and Businesses (ARIB). 

Direct-sequence (DS) spread-spectrum modulation is commonly used in CDMA 
systems, in which each information symbol is represented by a number of "chips". Representing 
one symbol by many chips gives rise to "spreading", as the latter typically requires more 
bandwidth to transmit. The sequence of chips is referred to as the spreading code or signature 
20 sequence. At a DS receiver, e.g., a rake receiver, the received signal is despread using a 

despreading code, which is typically the conjugate of the spreading code. IS-95 and J-STD-008 
are examples of DS CDMA standards. 

In the mobile radio channel, multi-path is created by reflection of the transmitted signal 
from obstacles in the environment, e.g., buildings, trees, cars, etc. In general, the mobile radio 
25 channel is a time varying multi-path channel due to the relative motion of the structures that 
create the multi-path. 

A characteristic of the multi-path channel is that each path through the channel may have 
a different phase. For example, if an ideal impulse is transmitted over a multi-path channel, 
each pulse of the received stream of pulses generally has a different phase from the other 
30 received pulses. This can result in signal fading. 

When multi-path propagation is present, the amplitude can vary dramatically. Multi- 
path propagation can also lead to time dispersion, which causes multiple, resolvable echoes of 
the signal to be received. In the receiver, correlators are aligned with the different echoes. Once 
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the despread values have been weighted, they are summed. This weighting and summing 
operation is commonly referred to as rake combining, 

FIG. IB illustrates a conventional radio receiver employing a channel estimator. The 
receiver includes an antenna 10 for receiving signals, a radio receiver 1 1 for filtering and 
5 amplifying signals and converting them to a suitable form for processing, such as complex 
numerical sample values, a channel estimator 12 which correlates received signal samples with 
known symbols stored or generated locally in a local code generator 14, to provide channel 
estimates, and a data decoder 13 for despreading and processing despread signal samples 
together with the channel estimates to extract information. Data decoder 13 may be, e.g., a rake 
10 receiver operating in a manner described in U.S. Patent No. 5,572,552 to Dent and Bottomley, 
which is hereby incorporated by reference. A channel estimate updater 15 updates channel 
estimates based on the latest data and/or pilot symbol decisions and the despread signal samples 
and provides this updated information to the decoder 13. Decoded data is output for further 
processing. 

15 Coherent detection requires estimation of how the signals were modified by the 

transmitter, channel, and/or radio processor. As discussed above, the transmission medium 
introduces phase and amplitude changes in the signal, as a result of multi-path propagation. The 
signal may also have become dispersed, giving rise to signal echoes, each echo having a phase 
and amplitude associated with it, represented by a complex channel coefficient. Each echo also 

20 has a delay associated with it. Coherent demodulation requires estimation of these delays and 
coefficients. Typically, the channel is modeled as discrete rays, with channel coefficients 
assigned to the different delays. 

Channel estimation for a received radio signal using both known modulation symbols 
embedded in the signal as well as unknown information symbols that are decoded by the 

25 receiver is described, for example, in U.S. Patent No. 5,335,250 to Dent et al., and also in U.S. 
Patents No. 5,331,666; No. 5,557,645; and No. 5,619,533 to Dent, all of which are hereby 
incorporated here by reference. Channel estimation specific to CDMA systems is described in 
U.S. Patents No. 5,151,919 and 5,218,619 to Dent, which are also hereby incorporated here by 
reference. 

30 More discussion of smoothing channel estimates using autoregression, that is IIR 

filtering, may be found in "A Wiener Filtering Approach to the Design of Tracking Algorithms 
with Applications in Mobile Radio Communications*', Ph.D. thesis of Lars Lindbom, Uppsala 
University (1995), which is also hereby incorporated by reference herein. This document 
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describes the benefit of adapting a smoothing filter's characteristics to the fading spectrum of the 
signal. 

In the past, the fading spectrum of a signal was assumed to be symmetrical. This is 
probably true in the long term (i.e., over several minutes), in accordance v^ith Jakes' model for 
5 fading in the urban, mobile radio propagation environment. More discussion of Jakes' model 
and modifications of the model to speed computation during simulations of communication 
system performance may be found in P. Dent, G. E. Bottomley, and T. Croft, "Jakes' Fading 
Model Revisited", Electronics Letters, vol. 29, no. 13, pp. 1162-1163 (June 24, 1993), which is 
hereby incorporated here by reference. 

10 Jakes' model assumes a uniform angular distribution of reflecting objects around a 

mobile receiver. The relative Doppler shift of reflected signals arriving at the mobile station at 
different angles relative to the direction of movement of the mobile station varies with the 
cosine of the angle of arrival. With a uniform angular distribution, the Doppler spectrum is then 
symmetrical and two sided, having as much reflected energy arriving ft-om behind the mobile 

1 5 station with a negative Doppler firequency shift as fi-om ahead of the mobile station, having a 
positive Doppler shift. Rays reaching the mobile station fi-om behind have clearly not 
propagated an equal distance from transmitter to receiver as rays reaching the mobile station 
from the front. However, these delay differences have typically been ignored, Jakes' model 
assuming that rays with such delay differences could nevertheless be combined to produce a net 

20 fading waveform for a path of delay equal to the mean of these rays. 

More specifically, delays lying within ±0.5 of a modulation symbol period of each other 
were combined to produce a net fading ray with a mean delay. Delays outside that ±0.5 
modulation symbol interval were grouped into a different ± 0.5 symbol window to obtain a 
different net fading waveform with a different mean delay. The different net fading waveforms 

25 with their associated modulation-symbol-spaced delays were then taken to characterize a 

multipath channel, each of the multiple paths nevertheless assuming to conform to Jakes' fading 
model, i.e., each path is the combination of rays arriving uniformly from all directions. 

In a CDMA system, particularly a wideband CDMA (WCDMA) system, chip intervals 
are much shorter, allowing multiple propagation paths to be resolved with much finer time 

30 resolution. Thus, it is no longer vaUd to use a Jakes' model which adds rays that differ in their 
propagation delays by even a fraction of a microsecond. This addition was valid only in the 
context of narrowband FDMA or medium bandwidth TDMA systems. In WCDMA systems, it 
is necessary to restrict the combination of different rays reaching the receiver to rays that have 
the same propagation delay from the base station to the mobile station, within ± 0.5 of a CDMA 
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chip duration. In the proposed IMT2000 system for next generation mobile telephony which is 
based on DS-CDMA, a frame has a duration of 10 milliseconds and is divided into 16 slots, 
each slot being divided into 2560 chips. Depending on the communication channel, 2560 chips 
are grouped into a number of symbols. For example, in the so-called Perch 1 Channel, there are 
5 ten symbols of 256 chips each. A certain number of these symbols are already knov^n and 
transmitted as pilots from the BS to MSs. One symbol in every paging channel slot is a so- 
called Perch 2 code. An exemplary CDMA signal format is shown in FIG. 2. 

In a 5 MHZ wide WCDMA system, a chip duration is typically 0.25 microsecond ()isec), 
so ± 0.5 chips has a duration of ± 0.125 ^sec, which may be expressed as ± 37.5 meters in terms 
10 of propagation path length variation. 

It may be shown that rays with the same delay to this order of accuracy must have 
reflected from objects lying on an elliptical contour having the base station and the mobile 
station as its foci. This may be understood with reference to FIG. 3 which depicts elliptical 
contours representing loci of objects from which the reflection delay is the same. For example, 
15 a ray reflected from objects lying on the elliptical contour closest to he mobile station have a 
delay of Tl, rays reflected from objects lying on the next closest elliptical contour have a delay 
of T2, rays reflected from object lying on the next elliptical contour have a delay of T3, etc. 

These objects are not uniformly distributed in angle around the mobile station, nor are 
they spaced at the same distance from either the mobile station or the base station. Moreover, 
20 since the base station lies inside the elliptical contour, if, as is usual, it employs directional 
transmit antennae, objects around the elliptical contour will not be uniformly illuminated. 
Consequently, the fading spectrum of a ray of given delay within ±0.5 chip periods will no 
longer be symmetrical about zero frequency. Rather, the fading spectrum of such a ray will be 
asymmetrical. This is illustrated in the power spectral plot shown in FIG. 4. 
25 In addition, the offset from zero frequency of the centroid of the fading spectrum is no 

longer independent of the direction of motion. 

One of the purported advantages of WCDMA is that the high time resolution enables 
resolution of individual reflecting objects such that each resolved ray is a single, non- fading ray, 
i.e., WCDMA is purported to eliminate fading. Of course, it is recognized that such "non- 
30 fading" rays will come and go, but on the relatively longer timescale of lognormal shadowing, 
which is easier to track. However, each ray will have a varying Doppler spread which means 
that its phase still varies at up to the Doppler rate, even if its amplitude varies much slower. 
Thus, there is still the need to track the varying complex value of the propagation channel in 
order to effect coherent signal decoding, i.e., with knowledge of a phase reference. Moreover, 
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the complete elimination of fading by resolving small reflecting objects is not achieved except 
using very large bandwidths, beyond the bandwidths of anticipated WCDMA systems, v^hich 
therefore find themselves in the intermediate region of propagation paths that still each comprise 
multiple rays. Fading models and channel estimation means for these WCDMA systems are 
5 addressed in U.S. Patent Application No. 09/227,180 filed January 7, 1999, in the name of Paul 
Dent, and entitled "Smoothing Channel Estimates by Spectral Estimation". This application is 
hereby incorporated here by reference. 

A receiver in w^hich the smoothing of channel estimates is adapted separately for per ray 
channel estimation according to the above-incorporated application is shovm in FIG. 5, and 

10 channel estimation may be performed as in FIG. 6. 

In a DS spread-spectrum receiver, a frequency offset or deviation may exist between the 
transmitter carrier frequency and the local oscillator of the receiver. The frequency offset results 
from different factors, including temperature variation, aging, and manufacturing tolerances. To 
address this offset, a phase ramp can be estimated and compensated for in an Automatic 

15 Frequency Control (AFC) control loop. Estimation can be based on a pilot channel, pilot 
symbols, or data symbols with a decision feedback. 

In the forward (base-to-remote) hnk of the IS-95 DS-CDMA system, a pilot channel is 
available for frequency offset estimation. The pilot is transmitted continuously, allowing 
tracking of variations in the offset. 

20 It is well known that a receiver AFC provides an estimate of the local crystal reference 

oscillator's error relative to the remote transmitter, and that this AFC estimate can be used to 
correct the crystal oscillator in order to correct the local transmitter frequency, as shown in 
FIG. 7. However, past attempts to ascribe observed variations in the received signal in part to 
Doppler channel variation via a channel estimate and in part to a crystal frequency error were 

25 inaccurate, as the channel estimate absorbed part of the frequency error. 

These problems are made even worse in a handover situation, in which a mobile station 
communicating with one base station at a certain frequency offset is handed over to another base 
station with another frequency offset. 

Therefore, there is a need to improve upon the discrimination of signal variations due to 

30 channel variations from signal variations due to oscillator error. 

SUMMARY 

It is therefore an object of the present invention to provide a technique for estimating an 
error between a carrier frequency of a transmitted signal and a local frequency reference in a 
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receiver that accounts for Doppler offset and crystal frequency error. It is yet another object of 
the present invention to provide a technique for estimating a frequency error between a local 
frequency reference of a receiver and carrier frequencies of multiple transmitters. 

According to exemplciry embodiments, this and other objects are met by a method and 
5 apparatus for processing code division multiple access signals received through a multipath 
propagation channel to control the frequency of a local frequency reference oscillator. The 
received signals are processed using the local frequency reference oscillator to obtain 
representative complex numerical samples for processing. The complex numerical samples are 
correlated with shifts of a locally generated despreading code. 

1 0 According to one embodiment, the correlated signals are complex channel estimates, 

each corresponding to a different ray of the multipath propagation channel. A frequency error 
estimate is computed for each ray based on successive values of a respective one of the channel 
estimates. A summer performs a weighted summation of the frequency error estimates to 
compute a relative frequency error estimate of the relative frequency error. The frequency error 

15 estimate may be computed by multiplying the current value of the respective channel estimate 
with the complex conjugate of a previous value of the same channel estimate and using the 
product as the frequency error estimate for the respective ray. An outer loop integrator may be 
used for integrating the frequency error estimates to produce a frequency error estimate and a 
control signal to control the local frequency reference oscillator to a value based on the received 

20 signal. 

According to another embodiment, the correlated signals are streams of complex 
despread values corresponding to each ray and successive symbol interval. Frequency errors on 
each of the despread value streams are corrected by progressively rotating the phase angle of 
successive despread values at a rate given by an associated frequency error integral. The 

25 frequency-corrected despread value streams are processed to produce complex channel estimates 
for each ray, and the residual frequency error estimate for each ray is determined by processing 
successive values of the channel estimates for the corresponding ray. The frequency error 
estimates are combined to produce a frequency error estimate and a control signal to control the 
local frequency reference oscillator to a value based on the received signal. Inner loop 

30 integrators may be used for integrating respective residual frequency error values to produce 
frequency error integrals, wherein the combiner combines the frequency error integrals to 
produce the control signal. Altemately, an inner loop integrator may be used for integrating the 
frequency error estimate to produce the inner loop integral values, and an outer loop integrator 
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may be used for integrating the inner loop integral values to produce a control signal to control 
the local frequency reference oscillator to a value based on the received signal. 

According to various aspects, a frequency error estimate may be computed for signals 
received from one transmitter or from multiple transmitters. For signals received from multiple 
5 transmitters, the frequency error estimates may be computed separately for each transmitter and 
then combined. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The features, objects, and advantages of this invention will become apparent by reading 
10 this description in conjunction with the accompanying drawings, in which hke reference 
numerals refer to like elements and in which: 

FIG. 1 A is a block diagram of an exemplary cellular radiotelephone communication 

system; 

FIG. IB illustrates a conventional radio receiver; 
1 5 FIG. 2 illustrates a format of a CDMA signal; 

FIG. 3 illustrates loci of reflecting objects of constant path delay; 
FIG. 4 illustrates a power spectral plot; 

FIG. 5 illustrates a receiver in which the smoothing of channel estimation is adapted 
separately for each ray; 
20 FIG. 6 illustrates a conventional method for channel estimation; 

FIG. 7 illustrates a conventional radio with transmit AFC from the receiver; 
FIGS. 8 A and 8B illustrate an apparatus and a method, respectively, for providing 
simultaneous frequency and channel estimates for each ray, according to one embodiment; 

FIGS. 9 A and 9B illustrate an apparatus and a method, respectively, for computing per- 
25 ray frequency error estimates by processing successive per-ray channel estimates according to a 
second embodiment; 

FIG. 10 illustrates an exemplary apparatus for providing a common frequency error 
estimate by combining per-ray despread values; 

FIG. 1 1 illustrates an exemplary apparatus for providing a common frequency error 
30 estimate using an AFC loop without an inner AFC loop; 

FIG. 12 illustrates an exemplary apparatus for providing a common frequency error 
estimate using an AFC loop with a common inner AFC loop; 

FIG. 13 illustrates an exemplary apparatus for providing a common frequency error 
estimate using an AFC loop with per-ray inner AFC loops; 
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FIGS. 14A-14C illustrate local frequency offsets in a handover scenario; and 
FIGS. 15A and 15B illustrate an exemplary apparatus and method, respectively, for 
estimating a frequency error in signals received from multiple transmitters. 

5 DETAILED DESCRIPTION 

For illustrative purposes, the following description is directed to a cellular radio 
communication system, but it will be understood that this invention is not so limited and applies 
to other types of communication systems. 

According to exemplary embodiments, methods and devices are provided for estimating 

10 a frequency error between a local frequency reference and at least one carrier frequency of at 
least one transmitted signal. 

According to one exemplary embodiment, the frequency error is computed by obtaining 
frequency error and channel estimates simultaneously. Referring to FIG. 8A, a received signal 
is despread in despreading units 300-1, 300-2, . . 300-n to obtain a number of despread values, 

15 corresponding to different rays. Despreading may be performed by correlating the complex 
received signal samples with time-shifts of a locally generated despreading code. An example 
of an efficient device for simultaneously correlating the same received signal with many time- 
shifts of a code is described in U.S. Patent No. 5,9631,893 to Dent and Wang, which is hereby 
incorporated here by reference. 

20 Despreading units 300-1, 300-2, . . ., 300 n also apply a per-ray frequency error 

correction, which is different for each ray in this illustrative embodiment. It will be appreciated, 
however, that the frequency error correction may be a common frequency error correction, i.e., 
the same correction may be applied for each ray, and may be applied once on the received signal 
input to blocks 300-1 to 300-n instead of separately at each block. When separate frequency 

25 corrections are applied, this results in a cumulative phase-untwisting of the phase angle of 
successive complex output values to compensate for the tendency of the frequency error to 
progressively rotate the phase angle. That is, if the frequency error on ray 2 that is despread by 
unit 300-2 causes a phase rotation of +10 degrees per output sample, then the frequency 
correction will rotate one despread output sample by -10 degrees, the next by -20 degrees, the 

30 next by -30 degrees and so forth. If the frequency correction is accurate, this eliminates the 

tendency of the output samples to progressively rotate in phase. To determine that this has been 
accomplished, per-ray frequency and channel estimators 12-1, 12-2, . . 12-n determine 
separately for each ray whether any progressive rotation tendency remains. Channel estimation 
may be carried out by observing the despread values when the underlying information symbol 
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modulation includes known symbols. For example, if the underlying information symbol 
modulation is QPSK, in which symbols have nominally unit amplitudes and phases of either 
±45 degrees or ±135 degrees, and a particular 2-bit symbol is known, then the despread value 
may be "unmodulated" by removing the phase angle modulation due to the known symbol, the 
5 remaining value being due only to the propagation channers phase and amplitude. By filtering 
or adding such values obtained from known symbols, an average value of the channel estimate 
for a particular multipath ray is obtained. 

In the signal format shown in FIG. 2, each of the slots of 2560 chips contains a number 
of known symbols. If each symbol is spread to 256 chips, then each slot contains 10 symbols. 

10 A certain number, e.g., four, of these symbols per slot may typically be known symbols used for 
channel estimation. Thus, channel estimates based on a certain number of these "pilot" symbols 
are obtained every slot, and their complex values should not exhibit a systematic phase rotation 
if frequency error has been removed. Consequently, the residual frequency error can be 
determined by measuring the difference in phase from one channel estimate to the next. This 

1 5 may be computed by forming the product: 

Cki . Ck* i. 1 

where Ckj is the channel estimate for ray k, slot i and Ck*-_i is the conjugate of the channel 
estimate for ray k in the previous slot. It will be recalled that channel estimates have had known 
symbol modulation removed. Therefore, the only cause of phase rotation from one slot to the 

20 next is uncorrected frequency error. This shows up as a non-zero imaginary part of the above 
product. The angle of this frequency estimate corresponds to a frequency error and can be used 
with appropriate smoothing to obtain an improved frequency error estimate. Although only the 
angle of the frequency estimate corresponds to the frequency error, in the interest of simplicity, 
the entire frequency estimate is referred to at various places in this application as the frequency 

25 error estimate. The magnitude of this differential phase estimate (or frequency estimate) is also 
proportional to the square of the ray magnitude | Ck| . In another implementation that will be 
described, the optimum common frequency error estimate across all rays is obtained by adding 
the imaginary parts of the above products across all rays (index k varying), as the squared- 
amplitude weighting thereby obtained is optimum for producing a common frequency error 

30 estimate. In FIG. 8A, however, the per-ray frequency error estimate from the imaginary part of 
the above product is applied to the ray-associated, AFC loop integrators 301-1, 301-2, . . 
301 -n, only one of which is shown for ease of illustration. The AFC loop integrator forms a first 
order AFC loop guaranteeing that the mean frequency errors from frequency error and channel 
estimators 12-1, 12-2 . . 12-n will be zero. If the mean frequency error is not zero, but, e.g.. 
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positive, the associated integrator ramps up and applies an ever increasing "detwisting" in the 
associated frequency error correction unit (one of 300-1, . . 300-n) until the mean frequency 
error is zero. At that point, the integrated frequency error from integrator 301-1 is the true 
mean frequency error for the associated ray. 
5 In addition to producing frequency error indicators, the device of FIG. 8 A includes a 

rake combiner 400 for producing soft data values. The rake combiner 400 multiplies frequency- 
corrected, despread values from despreaders 300-1 to 300-n with the conjugates of the 
corresponding channel estimates from channel and frequency error estimators 12-1, . . 12-n in 
complex multipliers 400-1, 400-2, . . ., 400-n. The results are then added in a summer 401 to 

10 produce a soft data value. The channel estimates are derived using known symbols, e.g., from a 
pilot channel, while the despread values relate to unknown symbols. Thus, the rake-combined 
soft values are indicative of the unknown symbols. The soft values from one whole frame of the 
signal format of FIG. 2 may constitute an error-correction and detection coded block. When the 
whole block is received, error correction and detection of the soft values is performed to obtain 

15 an error indication. If the error indication indicates correct decoding, then it is deemed that the 
received signal was a genuine signal and of sufficient reliability that the frequency errors present 
on the output of per-ray AFC loop integrators 301-1, 301-2, . . ., 301-n may be used to 
determine a correction to the crystal reference oscillator. This overall frequency error estimate 
is determined from the per-ray frequency error estimates using the combiner 402. 

20 While FIG. 8 A shows a Rake Combiner 400 for processing soft values, one skilled in the 

art will appreciate that other forms of demodulation may be used, including forms of multiuser 
detection, as well as the approach described in the co-pending U.S. Application No. 09/165,647, 
filed October 2, 1998, in the name of G. Bottomley. This application is hereby incorporated 
here by reference. 

25 There are many ways to determine an overall frequency error estimate, also referred to as 

a relative frequency error estimate, from the per-ray frequency errors, all of which may be 
suitable for implementing this invention. The techniques for determining the overall frequency 
error estimate differ in how much weight is placed on the per-ray frequency error estimates 
depending on their signal-to-noise ratios. In one example, no weighting is used. The rays 

30 processed in FIG. 8A have already been assumed to have been selected to be the strongest N 
rays out of total number of tested ray positions greater than N, and thus their frequency error 
estimates are simply averaged in the combiner 402. In another example, it can be argued that 
three positive frequency errors of +100 Hz, +110 Hz and +120 Hz and one negative frequency 
error of -130 Hz are just indicative of a Doppler spread between -130 Hz and +120 Hz, with 



-11- 



Attorney's Docket No. 040070-423 
Patent 

more reflecting objects giving rise to multipath rays of positive Doppler. According to this 
argument, the overall frequency error is half the sum of the extreme positive value of +120 Hz 
and the extreme negative value of -130 Hz, giving a result of -5 Hz. The third option is to 
weight the individual frequency errors by the average energy of the associated ray, which, for 
5 example, is available as the real part of the product Ck^ . Ck* j. 

In all cases, if the overall frequency error only needs to be estimated every 10 msec 
frame, the quantities from which it is derived should be averaged over that time period to reduce 
noise. Alternatively, the overall frequency error can be derived more often and averaged over a 
10 msec frame period to reduce variance. 

10 FIG. 8B illustrates an exemplary method for estimating a frequency error according to 

the first embodiment. The method begins at step 800 at which complex despread value streams 
are produced. At step 810, frequency errors for each stream are corrected. At step 820, channel 
estimates are produced. At step 830, the residual frequency error for each ray is estimated. At 
step 840, the frequency error estimates are combined to form a relative frequency error estimate. 

1 5 It may not be necessary in a particular application to use per-ray frequency correction. 

Usually, the channel estimators are sufficiently agile to track ray-to-ray and slot-to-slot phase 
variations, as long as the overall frequency error is not large. Thus, it is perhaps of greater 
importance to ensure a low overall frequency error. 

Thus, according to a second embodiment, a device is provided that uses per-ray channel 

20 estimation but not per-ray frequency correction. Referring to FIG. 9A, a device according to the 
second embodiment includes a per-ray channel estimator (12-1, 12-2, . . 12-n) and a rake 
combiner 400 as in FIG. 8A, but there is no per-ray frequency correction. Frequency error 
estimators 302-1, 302-2, . . 302 -n for each ray determine any systematic phase rotation 
tendency for each ray by multiplying a channel estimate at time (i) by the conjugate of the 

25 channel estimate at time (i-1). 

The sum of these products is then accumulated in the combiner 402, which is labeled the 
same in FIG. 9A as in FIG. 8A to indicate similarity of function. The sum from the combiner 
402 has an imaginary part and a real part, and the overall frequency error estimate is 
ARCTAN(imag/real) (or (imag/real) for small errors). The imaginary and the real parts may be 

30 averaged over the estimate period (e.g., one 10 msec frame) before computing their ratio or 
ARCTANGENT. As an alternative, frequency error estimators 302-1, 302-2, . . ., 302-n may 
produce a per-ray frequency error estimate, by multiplying a channel estimate at time (i) by the 
conjugate of the channel estimate at time (i-1) and then computing an ARCTAN (imag/real). In 
this case, combiner 402 does not require an ARCTAN operation. 
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The average result is then dumped to a control processor such as the control 
microprocessor 19 of FIG. 7, which, according to an exemplary embodiment, contains a 
software AFC loop integrator to form the overall AFC loop including the crystal reference 
oscillator. The software AFC loop integrator can, for example, comprise an integer variable to 
5 which the new frequency error estimate is added every 10 msec frame. The accumulated 

frequency error estimate is then used to output a control code via a digital-to-analog converter to 
obtain an analog control voltage for a voltage-controlled crystal oscillator, thereby adjusting its 
frequency until the mean of the successive frequency error estimates is zero, and the loop 
integrator reaches a more or less stable mean value. In a refinement, the stable value can be 

10 paired with a temperature estimate and stored by the control processor 19 in a crystal 

temperature-compensating table, from where it can be retrieved as an initial value upon next 
using the apparatus after a period of shut down. Entry of values into such a long-term 
temperature table can be severely restricted by various sanity checks and verification that the 
signal being received was from an authorized base station transmitter. 

15 In FIG. 9A, the combined or relative frequency correction is used to correct the crystal. 

In addition, or altematively, the combined frequency error correction may be used to rotate the 
channel estimates for use in rake combining, as fiirther discussed in co-pending U.S. 
Application No. 08/991,770, filed December 16, 1997, in the names of G. Bottomley et al. This 
application is hereby incorporated here by reference. 

20 FIG. 9B illustrates an exemplary method for estimating a frequency error according to 

the second embodiment. The method shown in FIG. 9B is similar to that shown in FIG. 8B, 
except that the frequency error estimate is computed at step 910 from the channel estimates 
produced at step 900, with no frequency correction taking place in advance. A weighted sum of 
the frequency error estimates is produced at step 920. 

25 In the devices shown in FIGS. 8 A and 9 A, the frequency error arid channel estimates are 

based on the use of pilot symbols embedded periodically in the transmitted signal format or on a 
common pilot channel. However, in some cases, an overall frequency error estimate should not 
be believed or acted upon unless the error correction and detection decoder indicates successful 
decoding of unknown symbols. Thus, with some increase in complexity, it is possible to defer 

30 making an overall frequency error estimate until the unknown symbols have been decoded and 
thus become known symbols. Substantially all symbols may then be used in a retrospective 
channel and frequency error estimation procedure which is otherwise identical to the process 
already described above. The "post-decoding" frequency error estimate will benefit from the 
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error correction decoding process and provide a more accurate and less noisy frequency error 
estimate. 

FIGS. 10-13 illustrate various implementations of an overall AFC loop according to 
exemplary embodiments. In FIG. 10, a common AFC loop including an AFC loop integrator 
5 508 develops an estimate of overall or relative frequency error. The frequency error from the 
integrator 508 is converted to a progressively rotating phase by successively adding the 
frequency error W to an accumulator 501, thereby generating Wt for progressively increasing 
time t. The accumulator 501 operates modulo-27t, i.e., upon overflov^, the remainder in the 
accumulator 501 continues to correctly represent the angle, now reduced by In, The cos/sin 

10 ROM 502 uses the instantaneous angle from the accumulator 501 as an address and outputs 

exp(-jWt) = cos(Wt) - jsin(Wt). This complex number is multiplied using respective multipliers 
503a- 1, 503b- 1, etc., by the despread outputs for each ray from the despreading unit 500 to 
apply the frequency correction prior to channel estimation in channel estimators 503a-3, 503b-3, 
. . ., 503n-3. Thus, the despread values are "untwisted" to remove systematic phase rotation 

15 caused by frequency error for the purposes of channel estimation. The untwisted despread 

values are also delayed by one information symbol period in per-ray delay units 503a-2, 503b-2, 
. . ., 503n-2, and the conjugates of the delayed values are multiplied by the corresponding 
undelayed values in multipliers 503a-4, 503b-4, . . ., 503n-4. The multiplied results are then 
summed in a summer 504, which yields a sum in which, if the effects of symbol-period to 

20 symbol-period changes in the information symbols were not present, the change in phase over a 
symbol period which represents uncorrected frequency error would appear as a non-zero 
imaginary part in the sum output. 

A unit 507 computes ARCTANGENT (imag/real) to yield the change in angle over the 
one symbol period of delay. This arctangent fiinction is a so-called two-argument arctangent 

25 which places the angle in the correct quadrant based on the signs of the real and imaginary parts 
out of summer 504, Then, the unit 507 subtracts the change of angle caused by the change of 
the modulating information symbol from symbol S(i-l) to symbol S(i). This of course 
presupposes that the symbol change is known. For example, any pilot symbols and/or unknown 
data symbols can be used to control frequency error by waiting until after symbols are error 

30 correction decoded and then retrospectively updating the AFC loop integrator or by making a 
tentative decision on the most likely value of the unknown symbol change based on the 
differential phase value. For example, if the differential phase value is within ±45 degrees, the 
differential change due to symbol modulation may be assumed to be zero S(i) = S(i-l), and the 
phase change may be ascribed to frequency error. On the other hand, if the differential phase 
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change is between +45 and +135 degrees, it may be assumed that the differential symbol change 
was responsible for 90 degrees of phase change. Then, 90 degrees is subtracted, and the residual 
differential phase change is used as the frequency error. The latter reduction of the phase 
modulo Tz/2 may be accomplished by choosing the appropriate modulo -In digital phase 
5 representation as a fixed point integer, adding 45 degrees, masking off the most significant two 
bits, then subtracting 45 degrees. 

The system of FIG. 10 is a self-contained receiver AFC loop that corrects a fi'equency 
error for channel estimation and data demodulation. However, it is also desirable to correct the 
local transmitter fi-equency by correcting the local crystal reference oscillator. The fi'equency 

10 error W present on the receiver AFC loop integrator 508 is thus output to the main control 
microprocessor 19 where it is digitally accumulated in another, slower loop integrator. 
Moreover, particularly when data symbols are used, the values are preferably accumulated in the 
digital (software) integrator only upon error correction and detection decoding indicating no 
uncorrected data errors are present. The accumulated value in the software loop integrator in the 

15 control processor 19 is then output via the D/A converter 17 of FIG. 7 to adjust a reference 

crystal oscillator such as the oscillator 16 of FIG. 7, thus closing a slower AFC loop around the 
reference oscillator 16. The closed loop transfer fimction of the inner (receiver) AFC loop 
including a loop integrator 508 acts like a low-pass pole within the outer AFC loop including the 
control processor 19, D/A convertor 17 and crystal oscillator 16. The design of the outer loop is 

20 stable and well-damped when this inner-loop pole is included. 

FIG. 1 1 illustrates the general form of a common AFC loop without per-ray AFC 
according to another embodiment. A duplexor 600 connects the antenna to both the receiver 
601 and the transmitter 602, which both use the crystal reference oscillator 16 for accurate 
firequency channel determination. The receiver 601 filters, down converts, samples and digitizes 

25 the received signal to obtain numerical samples at a multiple of the CDMA chip rate in a 

conventional manner. A despreader 603 uses local CDMA code generators in correlation with 
the numerical samples to despread the CDMA samples to obtain one complex number per 
information symbol for each of the selected multipath rays of significance. Unmodulators, one 
for each ray, remove phase changes due to known symbol modulation (if necessary). If the pilot 

30 symbols are all the same, viewing the common pilot channel as a series of pilot symbols, there is 
no need to remove the modulation. Unmodulation may be performed by multiplying the 
symbols by their conjugates. The unmodulated samples are accumulated or low-pass filtered 
over some period of known symbols. It will be appreciated that, when the known symbols, such 
as pilot symbols, are inserted at the remote transmitter, the combination of unmodulating and 
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filtering is equivalent to pilot-assisted, per-ray channel estimators 12-1, . . 12-n. However, 
data symbols that only become known after decoding may also be used for unmodulation. 

The filtered, unmodulated, per-ray complex values are then subject to differential phase 
discrimination (302-1, . . ., 302-n) which multiplies the latest value by the conjugate of a 
5 previous, delayed value in order to determine a complex number related to the phase change 
over the delay interval. The delay interval may be one symbol period, one slot period (referring 
to FIG. 2) or just one or more intervals between successive outputs of the channel estimator. 
However, the interval should not be so long that the phase change in the period could suffer 
from modulo-271 ambiguity or "wraparound". The choice of the interval therefore determines a 

10 maximum frequency error that can correctly be estimated. The outputs of the differential phase 
discriminators are already correctly weighted to be added. The sum of these differential 
complex numbers is formed in a complex summer 504. It will be appreciated that the single 
arrows entering and leaving boxes such as 302-1 or 504 or 507a are in fact carrying complex 
numbers that have both a real and an imaginary part. The complex number output from the 

15 summer 504 may be optionally ftirther low-pass filtered in a low-pass filter 507a, which 

includes a filtering channel for the real part and an identical filtering channel for the imaginary 
part. In the numerical domain, such filters can be finite impulse response (FIR) or infinite 
impulse response (IIR) filters defined by a number of coefficients that set the filter bandwidth or 
time-constants. 

20 The filtered, combined differential complex values are then converted to a phase change 

value over the delay interval, e.g., using two-argument ARCTAN function (i.e., four-quadrant 
ARCTAN) 507b. The phase change in a fixed delay interval is an estimate of the relative 
frequency error and is conveyed to a control microprocessor 19 where it is integrated in a 
software loop integrator 508, depending, e.g., on whether the microprocessor receives a separate 

25 confirmation of correct decoding of received data from an error correction/detection decoder. 
The latter criterion is not intended to be a limitation of the invention but represents a feature of 
the preferred implementation. The software loop integrator value is then output to the DtoA 
convertor 17 to form an analog control voltage for the crystal reference oscillator 16, thereby 
correcting frequency error for both the receiver 601 and the transmitter 602. The use of a D/A 

30 converter 17 driving a voltage-controlled crystal oscillator (VCXO) is not meant to be a 

limitation of the invention, but is a typical practical implementation. Alternately, a frequency 
synthesizer with fine frequency control steps, such as a Direct Digital Synthesizer (DDS) or a 
fractional-N synthesizer, may be used, in which case the software loop-integrator 508 output 
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value would be added to the synthesizer control code to effect fine adjustment of the receiver 
601 or transmitter 602 channel frequency. 

FIG. 12 shows yet another overall AFC loop including an inner AFC loop to adjust the 
receiver fi-equency error and an outer AFC loop to correct the transmitter frequency based on the 
5 correction developed by the inner loop for the receiver. For brevity, only the differences from 
FIG. 10 will be highlighted. The inner AFC loop comprises using a numerically controlled 
oscillator (NCO) 720 that untwists a progressively rotating signal phase, which can be applied 
by complex multiplying the signal either before or after despreading in a despreader 603 with 
the NCO's output of exp(-jWt). The NCO 603 performs the same function as the modulo-27t 

10 accumulator 501 and cos/sin ROM 502 of FIG. 10. The fi-equency control value to NCO 720 is 
developed by inner loop integrator 510 which integrates the combined firequency error estimates 
fi-om ARCTAN circuit 507b. This inner loop operates until the output of inner loop integrator 
510 is equal to the correct mean or relative frequency error between the received signal and the 
receiver. This frequency error is then applied to the outer loop integrator 508, which may be 

15 part of the control microprocessor 19. Differences between the inner loop and the outer loop 
can be, e.g., that the inner loop operates more rapidly to correct receive frequency error without 
waiting for the error correction decoder to declare error free data decoding. The outer Joop 
operates more slowly to provide long term correction of crystal oscillator 16 and is preferably 
protected from operating on erroneous data. 

20 FIG. 13 is similar to FIG. 12, except that the inner loops are separate for each multipath 

ray. In this embodiment, there is one NCO per ray, NCO 720-1, . . ., NCO 720-n. Likewise, 
there are per-ray inner loop integrators (510-1, .. ., 510-n) and a per-ray ARCTAN circuit, 
which is now generalized to Cartesian-to-polar converter 507b-l, . , ., 507b-n. The Cartesian-to- 
polar converter produces both an angle output and a signal strength output. The signal strength 

25 output can be related to the amplitude of the differential phase discriminator output, which is in 
turn related to the power of a ray. The inner loop thus now comprises a collection of n loops, 
one per significant multipath ray, and the respective loop integrators settle to output mean values 
related to the mean frequency error of the associated ray. The combiner 507 performs a 
weighted sum of the inner loop integrator outputs using the amplitude-related output of the 

30 associated Cartesian-to-polar convertor as the associated weight. The weighted output is then a 
estimate of the common or relative frequency error across all rays and is fed to the outer AFC 
loop integrator 508 included in control processor 19. The sum of the signal strength output of 
Cartesian-to-polar convenors 507b-l, . . ., 507-n may also be formed by a second summer (not 
shown) and used for automatic gain control to keep the sum constant at a desired value, for 
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example unity, in which case it is unnecessary to divide the combined output of 507 by the sum 
of the AFC weights. 

Various inventive configurations have been shown above for correcting the reference 
frequency oscillator of a transmitter-receiver based on fr'equency estimates made on separate 
5 multipath rays, as are usually detected in a receiver for CDMA signals. The various 

configurations can include a single, common AFC loop to correct the reference oscillator, a 
separate inner AFC loop to correct received signal frequency error, which correction is then 
applied to a slower outer loop to correct the reference crystal oscillator, separate inner AFC 
loops for each multipath ray, the per-ray inner-loop corrections then being optimally combined 

10 to form a common correction to the reference oscillator, or both per-ray and a common receiver 
AFC loop, the common receiver AFC loop then being an outer loop to the per ray loops. An 
outer loop can be added to further integrate the common receiver AFC loop signal to form a 
crystal reference oscillator correction loop for the transmitter. The latter configuration 
represents loops within a loop within a loop, and such configurations are allowed so long as 

1 5 every outer loop is properly designed to tolerate the closed loop poles of the inner loops within 
it while maintaining desirable stability and damping. 

The embodiments described above address how to handle a frequency error between a 
local frequency reference of a receiver such as an MS and the carrier frequency of a single 
transmitter, such as a BS. The frequency error problem is aggravated if the receiver receives 

20 signals from multiple transmitters at the same time, e.g., if a mobile station is communicating 
with more than one base station simultaneously, as in a soft handover situation. 

According to the WCDMA specification, two base stations are allowed to be frequency 
misadjusted by ±0. 1 ppm which corresponds to ±200 Hz at a carrier frequency of 2 GHz. At the 
same time, the MS can have frequency inaccuracy from the BS by ±0.1 ppm which also 

25 corresponds to ±200 Hz. This is described in more detail in "UTRA (BS) FDD; Radio 
Transmission and Reception", Technical Specification No. 3G TS 25.104, ver. 3.1.0, 3G 
Partnership Project (Dec. 1999), which is hereby incorporated here by reference. 

Suppose that the MS has a frequency offset equal to 0. 1 ppm and is going to camp to 
new BS which has a frequency offset of -0.1 ppm to the first BS. If the MS combines the 

30 received path-rays as described above to estimate the frequency offset, it is possible that the MS 
may not be able to adjust its local frequency to the new BS and thus may begin transmitting to 
this BS with a frequency offset that is more than ±0.1 ppm. The likelihood of this occurring 
increases if the number of path-rays received from the new BS is less than that received from the 
oldBS. 
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FIGS. 14A-14C illustrate the problem. As shown in FIG. 14A, the MS local frequency 
can be up to 200 Hz away from the BSl. In soft handover, when a new BS (BS2) with a 
frequency offset of -200 Hz from Bl, starts to communicate with the MS, it has a frequency 
offset of 400 Hz to the MS. As shown in FIG. 14B, if the MS employs the received path-rays 
5 from BSl and BS2 to estimated the frequency offset, the MS will estimate a frequency offset 
that is more than 200 Hz from the local frequency of BS2. Hence, the MS should adjust its local 
frequency to be in between the frequencies of BSl and BS2, in order to have a frequency offset 
in between ±200 Hz. 

This problem is addressed, according to another aspect of the invention, by handling 

10 received path rays from different BSs separately. According to an exemplary embodiment, the 
received path-rays from different BSs are labeled. The frequency offset is then estimated with 
respect to both BSs, despite the number of received path-rays from them. 

To keep the rays from different base stations separate, the received path rays from 
different BSs are labeled. The labels identify from which BS the received path-rays are 

15 transmitted. For example, in a WCDMA system. Long Codes may be used to identify the base 
stations from which signals are transmitted. The Long Codes may be applied and detected in 
any known manner, e.g., as described in U.S. Application No. 09/1 12,689, filed July 9, 1998, 
and hereby incorporated here by reference. 

The AFC algorithms described above can individually be applied to each BS, and the 

20 frequency offset between the MS and the BSs can be determined. A decision can then be made 
as to what the final frequency offset is by using those estimations so that the MS local frequency 
is in between the references from BSl and BS2, as shown in FIG. 14C. 

According to this other aspect, the MS can handle multiple BSs by using a separate inner 
loop AFC for each ray, irrespective of the BS from which the ray originates. Altematively, in 

25 FIG. 15A, received rays are grouped into a first group of rays originating from a first BS and a 
second group of rays originating from a second BS, based, e.g., on the Long Codes. Although, 
in the interest of simplicity, the device shown in FIG. 15A applies to two base stations, it will be 
appreciated that the invention is not limited handling signals received from two base stations but 
may be applicable to signals received from any number of base stations. 

30 Channel estimation, frequency error estimation and combining of frequency error 

estimates may be performed on the rays, e.g., as described in any of the embodiments above. 
The combined frequency error estimates in each group of rays are integrated in inner loop AFC 
integrators 508-1 and 508-2 and used to correct received frequency errors for the respective 
group of rays. Then the integrated frequency errors for all base stations are combined in 
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summer 517 to produce a relative frequency error estimate. The summed output is then 
integrated using an outer loop integrator 508 which may be a software AFC integrator. The 
outer loop integral is then applied to control the crystal oscillator. 

The BS-related frequency error integrals may be combined in a number of different 
5 ways, with or without weighting. For example, the combining in 517 may be an unweighted 
average, or each of the combined estimates may be weighted by the received signal strength 
from the respective BS. Altematively, the estimate from the BS that exhibits the strongest 
signal strength may be weighted 1, and the estimates from the BSs that exhibit the weaker signal 
strengths may be weighted 0. According to another combining method, only the inner loop 

10 AFC values associated with rays of the BS to which the mobile is transmitting, as determined by 
an assigned transmit spreading code, are fed to outer loop integrators. It is also permissible to 
ignore whether a ray originates from a particular BS and use the previously described combining 
methods, regardless of where a ray originates. 

FIG. 15B illustrates a method for estimating a relative frequency offset according to the 

15 second aspect. The process begins at step 1500 at which a determination is made which 
transmitter the received signal is transmitted from, based, e.g., on the Long Codes used for 
correlation. At step 1510, the frequency errors are estimated separately for each transmitter. At 
step 1 520, the frequency error estimates are combined. 

It will be appreciated by those of ordinary skill in the art that this invention can be 

20 embodied in other specific forms without departing from its essential character. The 

embodiments described above should therefore be considered in all respects to be illustrative 
and not restrictive. For example, although described above with reference to a CDMA 
communication system, the invention is also applicable in other types of communication 
systems. 
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